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We use the Lanczos method to calculate the variance a^(E,4>) of the number of energy 
levels in an energy window of width E below the Fermi energy for non-interacting dis- 
ordered electrons on a thin three-dimensional ring threaded by an Aharonov-Bohm flux 
<p. We confirm numerically that for small E the flux- dependent part of a^(E,<l>) is well 
described by the Altshuler-Shklovskii-diagram involving two Cooperons. However, in the 
absence of electron-electron interactions this result cannot be extrapolated to energies 
E where the energy-dependence of the average density of states becomes significant. We 
discuss consequences for persistent currents and argue that for the calculation of the 
difference between the canonical- and grand canonical current it is crucial to take the 
electron-electron interaction into account. 

The mesoscopic persistent current in a small metal ring threaded by an Aharonov- 
Bohm fhix (/) has been predicted long tirne ago by Hundlil. However, the role of 
disordem and electron-electron intoractionsS has only been addressed (a uite recently, 
and is not completely understoodaa'El. Jihc experiment by Levy et alB and the sub- 
sequent experiments by other groupsatl have motivated many recent theoretical 
works. There exists now general agreement that the surprisingly large experimen- 
tally observed persistent currentslH can only be explained by taking the electron- 
electron interaction into account. Nevertheless, the calculations within models of 
non-interacting electrons have lead to new insights into the nature of mesoscopic 
disordered systems. In particular, in a seminal paper Altshuler, Gefen and Imrjild 
have shown that for non-interacting electrons in the diffusive regime the dominant 
contribution to the average persistent current arises from striking differences_be- 
tween the canonical and grand-canonical ensembles. Because experimentallyQ no 
external leads were attached to the rings, one should calculate the canonical persis- 
tent current 

where F(N, ip) is the canonical free energy, the flux tp — (/>/0o is measured in units of 
the flux quantum (jjQ = hc/e, and — e is the charge of the electron. Unfortunately, the 
usual methods of quantum statistical mechanics are based on the gtand canonical 
formalism, where the chemical potential /i is held constant. In RefJlj the disorder 
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averaged canonical current I{N, ip) was therefore related to the variance 
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^'{^i\^)^N^^^*,^)-[Ni^^*,^)\ (2) 

of the particle number N(fj,, ip) in a corresponding grand-canonical ensemble at a 
certain value ^* of the chemical potential. Here and below the overbar denotes 
averaging over the disorder, fi* should be chosen such that the disorder- and flux- 
averaged particle number in the corresponding grand-canonical ensemble agrees 
with the given particle number N in the original canonical ensemble, i.e. 



dipN{^l*,ip)^N . (3) 







The expression derived in Ref.li can be written asliiJ 



where /gc (/i* , </?) is the average grand-canonical current, and K{fi, if) = {ON {jjt, Lp) / d 
is the average compressibility. For non-interacting electrons in the diffusive regime 
Igc{l^*Tf) is exponentially small; the leading interaction contributioiLto /gc(/^*,V') 
is not negligible and has been calculated by Ambegaokar and Eckerrfl. 

Note that Eqjl defines ^* as function of N , so that the right-hand side of Eq.^ 
is indeed a function of N and p. As it stands, Eq.^ is valid even for interacting 
systems, and can be used to take the effect of electron-electron intera£ti(ins on the 
canonical current in a simple but non-perturbative way into accountli3ilj. In this 
work we would like to restrict ourselves to non-interacting electrons in the diffusive 
regime. Then K~^(^*,Lp) can be approximated by a flux-independent constant A, 
the average level spacing at the Fermi energy /i*. In this way the calculation of 
the average canonical persistent current is reduced to the problem of calculating 
the flux- dependent part of the variance ip) of the particle number in a grand 

canonical ensemble. To calculate this quantity, let us write 

S^(M,'y5) = lim cr^(£',^,(^) , (5) 

with 

a''{E,ti,p)^ r de r de'K^ie.e'.if) , (6) 



K2{€,£',ip) = p{e,'fi)p{e',ip) - p{e,ip) p{e\ip) . (7) 

Here p(e, ip) = 5{e — ea{p)) is the density of states for a given realization of 
the disorder. (The £0(1^9) are the exact solutions of the Schrodinger equation for 
fixed random potential.) Note that cr^(i?, /i, (p) is the variance of the number of 
energy levels in an interval of width E below the Fermi energy. Combining Eqs.^-0, 
we obtain for the average canonical persistent current of non-interacting electrons 
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(ignoring the contribution /gc (/i* , (p) , which is cxponcntiahy small in the diffusive 
regime) 

-eA 9 



/(iV,^)«— -— hm / de de'K2{e,e',^) . (8) 

h 2 dip E^ccJ^,^j^ Jf,'-E 



This expression has been used in Rcfs.EilEJIl^ to calculate the average persistent 
current. For the covariance function 7^2 (e, e', (^) these authqis have substituted 
a perturbative approximation due to Altshuler and Shklovskiillj, which for a thin 
quasi one-dimensional ring (where diffusion is only possible along the circumference) 
is given by 
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27r2 ^ 1 [47r2i;c(fc + 2v9)2-i(e-e')+r]2 
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[Air^Eck^ -i{e-e') + TY 



(9) 



The first term in the curly braces is due to thii two-Cooperon diagram, while the 
second term is due to the two-Diffuson diagramllB. Here Ec is the Thouless energ5it3, 
and r is a cutoff energy that has been introduced by hand into the non-interacting 
model. For non-interacting electrons non-perturbative effects due to higher-pri 
terms give rise to contributions to F of the order of the average level spacinglJ 
In the presence of electron-electron interactions, F should take inelastic processes 
approximately into accounto. Eq.^ is believed to be the dominant contribution to 
K2 in the regime A;S|e — e'|^E?i/T, where r is the elastic lifetime. 

It is tempting to substitute Eqjp into Eq.^ and obtain in this way a simple an- 
alytic result for the average canonical persistent current00'El. However, such a 
procedure is based on a hidden assumption, which apparently has not been noticed 
in the recent literature on persistent currents: In the derivation of Eq.^ it is im- 
plicitly assumed that in the energy window of interest the energy- dependence of the 
average density of states can be neglected. On the other hand, according to Eq.|| 
the limit E ~* 00 has to be taken in order to obtain the physical persistent current, 
so that we need to know the function K2{e,e',ip) for all energies e and e' below 
the Fermi energy. Obviously, in this regime the energy-dependence of the average 
density of states p(e, if) cannot be ignored. Of course, if we proceed by performing 



the above substitution anyway, it is easy to show that (see Eq, 14 below) 



— if) ^—a^E,ti,p>) , foTa,\lE>E, . (10) 

Thus, at the first sight it seems that the limit i? 00 in Eq.^ is trivial and the 
above substitution is justified. However, this way of reasoning is clearly not self- 
consistent, because the result |l^ depends via Eq|9| on the assumption of a constant 
density of states. In this work we shall carefully examine this point numerically. 
Our main result is that for non-interacting electrons Eq.hd is not correct, so that 
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the physical persistent current cannot be calculated by simply substituting the 
Altshuler-Shklovskii result ^ into Eq^. We also present a numerical test of Eq.||. 
In a three-dimensional system Altshuler-Shklovskii scaling has recently been repro- 
duced numerically by Braun and MontambauxO, but these authors examined only 
the flux-independent part of (T^(£^,/i, ip). 

Using the Lanczos method, we have numerically calculated the exact energv 
levels Ea {<p) of the spinless nearest- neighbor tight-binding Anderson HamiltonianEd 
with diagonal disorder. The site-diagonal random potentials are assumed to be 
box-distributed in the interval [—w/2,w/2]. The Aharonov-Bohm flux is taken 
into account by choosing flux-dependent hopping energiesEil in the ±x-directions, 
t±x — te^'^'^^v/N^ ^ where is the number of lattice sites in the direction x. For 
convenience, all energies will be measured in units of the hopping energy i.e. we 
formally set t = 1. For a quasi one-dimensional ring with — 20, Ny — Nz — 5 
we found diffusive behavior for w = 2.5. However, in order to reduce the statistical 
errors so that even the flux-dependent part of a^{E, ^,ip) can be resolved, it is 
necessary to average over an extremely large number (typically 20000) of statistically 
independent realizations of the random potential. Our calculations were performed 
with the help of a parallel codcEl on a cluster of up to 30 work-stations. 

Because the Altshuler-Shklovskii result ^ is based on the assumption that in the 
energy window of interest the average density of states p(e, ip) can be approximated 
by a constant, it is important to choose the chemical potential /i properly. In Fig.^ 
we show the numerical result for the average density of states for a typical set of 
parameters. Obviously fi — —0.8 is a good choice, because then the Fermi energy 




Fig. 1. Average density of states for a 20 X 5 X 5 lattice. We used w = 2.5 and averaged over 
20000 realizations of the disorder. Thick line: (p = 0; thin line: (p = 1/4. The inset shows part of 
the central regime on a larger scale, together with our choice of the chemical potential fj, = —0.8. 
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lies on the right-hand side of a broad maximum, so that the width of the interval 
[jj, — E, /i] where p(e, tp) can be approximated by a constant is as large as possible. 

Next, let us test the Altshuler-Shklovskii prediction for the variance a'^{E,iJ,, (p) 
of the number of energy levels in an energy window of width E below the Fermi 
energy. For better comparison with the numerical results, it is convenient to expand 
the theoretical prediction for a-^{E, ii,(p) into a Fourier series. Substituting Eq.^ 
into Eq.0, we obtain 



a'^{E, n, ip) = al{E) + ^ 'jlmiE) cos(47rm(p) 



(11) 



where for m > 1 



<n{E) 



and the flux average is 



exp 
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X cos 
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m 
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, (12) 



(13) 



Here E = E/Ec, and T = T/Ec- Note that Ec (and hence E) implicitly depends on 
fj.. For £; > 1 > f Eqs.|l| and P reduce to 



■ exp 



4iE) ^ 4e'^' 



4ln[f-V2] 



1,2,. 



(14) 
(15) 



Hence, the higher harmonics cr|„j(£') become independent of E as soon as iJ 1, 
implying the validity of Eqjl^. Physically this means that the dominant contri- 
bution to the flux-dependent part of a'^{E, ^,ip) is due to states with energies in 
an interval of width Ec below the Fermi level. In Fig.^ we show our numerical 
results for the flux-dependent part of cr^{E, fj,, ip) for different values of i? = E/Ec- 
The flux-oscillations with a fundamental period of (/'o/2 are clearly visible, and the 
amplitude is in perfect agreement with the diagrannnatic calculation based on the 
two-Cooperon diagram of Altshuler and Shklovskiilij. Note also that the amplitude 
of the flux modulation does not change when the energy E is increased beyond 
the Thouless energy. To the best of our knowledge, this is the flrst numerical 
confirmation that the flux-dependence of cr^LE, f) is indeed correctly described 
by the two-Cooperon diagram given in Ref.E3. Our numerical data for the flux- 
independent part a'^{E, ij,ip) is shown in Fig.^. Although the overall quantitative 
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Fig. 2. The variance a^^ (E, fi, ip) as function of ip for different values of E = E/Ec for /i = —0.8 
and Ec = 0.01. From bottom to top: E = 1; E = 5; E = 15; E = 60. The solid lines are 
theoretical curves according to Eq.hd with F = 0.55. We have chosen the flux-average <Jq{E) as a 
fitting parameter. See Fig. for a numerical calculation of o"g(i?). 
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agreement is not as good as in the case of the flux-dependent part, the quasi-one- 
dimensional -n/E- behavior is clearly visible. It should be kept in mind, however, 
that the Altshuler-Shklovskii calculation is only valid for energies E ^ Ti/t. Using 
simple second order perturbation theory, we estimate H/t fa 0.5 in the parameter 
regime relevant to the figures. With Ec ~ O.OI, the theoretical curve in Fig.^ should 
be quantitatively accurate for i? ^ 50. 

The theoretical results |l2[-|l5| have been derived under the assumption that in 
the energy interval of interest the average density of states can be approximated by 
a constant. The important question is now whether it is allowed to use the small- -B 
regime as a basis for the extrapolation of the flux- dependent part of cr'^{E, /z, (p) to 
E ^ oo. According to Eq.|| this_extrapolation is necessary to obtain the persistent 
current. Note that in Refs]i3^llSlla such a procedure is adopted without further 
comment. In Fig.^ we have plotted our numerical result for the second harmonic 
(T|(i?) as function of energy. Although for sufficiently small E we obtain excehent 
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Fig. 4. Second Fourier component (Tj of '^'^{E, fJ., ifi) as function of E. The parameters are chosen 
as in Figs.0 and ^. The solid line is the theoretical prediction |l2| with f = 0.55. The upper labels 
on the horizontal axis and the right vertical scale refer to the corresponding average density of 
states p(e) at ip = (dashed line). 



agreement with the theoretical result at £' w 150 the numerical data begin 
to deviate strongly from the theoretical curve. Note that this energy cannot be 
identified with h/{TEc) ~ 50 - instead, it is the energy scale where the average 
density of states (shown as dashed line in Fig.^ starts to deviate significantly from 
its value at the Fermi energy. In contrast, according to the theoretical prediction 
the non-zero harmonics cr|^(i?) should saturate for E ^ 1, and nothing special 
should happen at i? = fi/ijEc) ~ 50. Indeed, at i? w 50 our numerical result 
for cr|(£^) still agrees well with the theoretical prediction. However, as soon as the 



8 Spectral correlations in disorderd mesoscopic metals and their relevance for persistent currents 



average density of states deviates from its value at the Fermi energy, our numerical 
data in Fig^ clearly disagree with the theoretical curve. We therefore conclude 
that the extrapolation of Eq.|l^ to ii^ ^ oo is not possible. Hence, at least for non- 
interacting electrons, the average persistent current cannot be calculated by simply 
substituting the Altshuler-Shklovskii result for a^{E, ii,ip) into the formula^ and 
taking the limit E ^ oo to obtain the variance ^p) of the total particle number. 

The physical reason is that S^(^, (^) depends on all energiesea (1^9) < ^. Because 
in the diffusive regime energy levels tend to repel each otheiO, the fluctuations of 
energy levels £0(1^9) deep inside the Fermi sea induce also changes in the energy 
levels in the vicinity of the Fermi energy. Thus, the variance of the total particle 
number depends on the statistical properties of the entire spectrum below the Fermi 
energy, so that the energy-dependence of the density of states cannot be ignored. 

Realistic models of non-interacting electrons certainly should have energy-dependent 
densities of states, so that for non-interacting electrons the flux-dependent part of 
the variance 5]^(/i,(/3) of the total particle number cannot be calculated from the 
Altshuler-Shklovskii diagram. In other words, in the absence of electron-electron 
interactions the second term in EqJJ is not determined by the statistical properties 
of the spectrum in the vicinity of the Fermi energy. On the other hand, in the 
presence of electron-electron interactions we expect a completely different scenario: 
in this case Eq.|l^ should be at least qualitatively correct, because in an interact- 
ing many-body system weakly damped propagating quasi-particles exist only in the 
vicinity of the Fermi energy. But only quasi-particles that can propagate coherently 
around the ring can probe the sensitivity to twists in the boundary conditions (and 
thus contribute to the flux-dependence of the particle number variance), so that 
we expect that electron-electron interactions will eliminate the contribution from 
many-body states with energies deep inside the Fermi sea. In other words the differ- 
ence I{N, (f) — Igc {fi* , (p) in Eq.^ should be extremenly sensitive to electron-electron 
interactions in the sense that only if the damping of the quasi-particles is taken into 
account, the flux-dependent part of S'^(/i*,(^) is determined by the spectrum close 
to the Fermi energy (see Eq.|o|). In this sense the procedure adopted in Refs.E^lliilll 
was physically correc although within a model of non-interacting electrons it 
cannot be justified. Note, however, in the presence of electron-electron interactions 
the contribution Jgc(/z*,v5) in Eq.^is not negligible, and it is not clear whether the 
persistent current is still dominated by the difference between the canonical and the 
grand-canonical current. 

We would like to emphasize that at this point theL_above scenario should be 
cons idered as a plausible hypothesis, which implicitlycJ has also been made in 
Refs.Ej'oli3. We suspect that in an interacting many-body system the microscopic 
mechanism which eliminates contributions from states deep inside the Fermi sea 
is closely related to the phenomenological cutoff parameter F in Eqj^. A pertur- 



bative calculation of 
reported by BlanterE 



in an interacting mesoscopic conductor has recently been 
but in the present context more accurate calculations are 



needed, which take also the energy-dependence of F into account. Although we have 
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started from a model of non-interacting electrons, we have arrived at the problem of 
electron-electron interactions. We believe that a satisfactory and generally accepted 
solution of the long-standing persistent current problem will only be obtained if the 
role of electron-electron interactions in mesoscopic disordered conductors is more 
thoroughly understood. 

To conclude, let us briefly summarize our main results: First of all, we have 
numerically confirmed that for a thin three-dimensional ring in the diffusive regime 
the variance cr^(i?, /i, (p) of the number of energy levels in a small interval of width 
E can indeed be described by the perturbative expression derived by Altshuler and 
Shklovskii0. Our second important result is the observation that for non-interacting 
electrons the flux-derivative ip)/dip of the variance of the total mimher of en- 

ergy levels below the Fermi energy cannot be obtained from the Altshuler-Shklovskii 
expression for a^{E, n, (p). We have argued that electron-electron interaction should 
modify this result, but a microscopic proof of this conjecture requires a better un- 
derstanding of dephasing in disordered mesoscopic conductors. We hope to come 
back to this point in a future publication. 

Acknowledgements 

This work was financially supported by the Deutsche Forschungsgemeinschaft (SFB 
345). We would like to thank J. A. Holm for providing us with his "Master-Slave 
Implementation for Parallel Virtual Machines" , which was crucial to run our nu- 
merical simulation efficiently on a work-station cluster. 

References 

1. F. Hund, Ann. Phys. (Leipzig) 32, 102 (1938). Prof. Hund is Professor Emeritus at the 
Institut fiir Theoretische Physik der Universitat Gottingen, and has recently celebrated 
his 100'^ birthday. 

2. M. Biittiker, Y. Imry, and R. Laudauer, Phys. Lett. 96 A, 365 (1983); 

3. V. Ambegaokar and U. Eckern, Phys. Rev. Lett. 65, 381 (1990); ibid. 67, 3192 (1991). 

4. P. Kopietz, Phys. Rev. Lett. 70, 3123 (1993); ibid. 71, 306 (1993) (E); P. Kopietz, Int. 
J. Mod. Phys. B 8, 2593 (1994). 

5. G. Vignale, Phys. Rev. Lett. 72, 433 (1994); P. Kopietz, Phys. Rev. Lett. 72, 434 
(1994). 

6. A. Ahland and Y. Gefen, Phys. Rev. Lett. 72, 2973 (1994); P. Kopietz, Phys. Rev. 
Lett. 72, 2974 (1994). 

7. L. P. Levy et ai, Phys. Rev. Lett. 64, 2074 (1990). 

8. V. Chandrasekhar et ai, Phys. Rev. Lett. 67, 3578 (1991). 

9. D. Mailly, C. Chapelier and A. Benoit, Phys. Rev. Lett. 70, 2020 (1993). 

10. B. L. Altshuler, Y. Gefen, and Y. Imry, Phys. Rev. Lett. 66, 88 (1991). 

11. There has been some confusion in the literature about the correct sign in Eq.^. Appar- 
ently this is due to the fact that in the present case we are dealing with three thermo- 
dynamic variables N, /i and (p which satisfy one relation of the form N — N{iJ,, p). 
Then we have (§^)Af/(^^)/i — "^^^ minus sign is familiar from elemen- 
tary thermodynamics, see E. M. Lifschitz and L. P. Pitajewski, Statistical Physics, 
Part I, (Pergamon Press, Oxford 1984), §16. If we approximate {^)n ~ and 



10 Spectral correlations in disorderd mesoscopic metals and their relevance for persistent currents 



(^)m ~ In ^^"^ cancel the 5(p, we obtain SN « —{^)^p6^. The minus-sign is 
missed if one naively treats 5N and as conventional differentials. 

12. R. Berkovits and Y. Avishai, Phys. Rev. Lett. 76, 291 (1996). 

13. A. Volker and P. Kopietz, preprint ;ond-mat/9612115| , accepted for publication 
in Z. Phys. B. 

14. A. Schmid, Phys. Rev. Lett. 66, 80 (1991). 

15. F. von Oppen and E. K. Riedel, Phys. Rev. Lett. 66, 84 (1991). 

16. B. L. Ahshuler and B. L Shklovskii, Zh. Eksp. Teor. Fiz. 91, 220 (1986) [ Sov. 
Phys. JETP 64, 127 (1986)]. 

17. K. B. Efetov, Adv. Phys. 32, 53 (1983). 

18. N. Dupuis and G. Montambaux, Phys. Rev. B 43, 14390 (1991); J. Fricke and 
P. Kopietz, Phys. Rev. B 52, 2728 (1995). 

19. D. Braun and G. Montambaux, Phys. Rev. B 52, 13906 (1995). 

20. P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

21. H. Bouchiat and G. Montambaux, J. Phys. (Paris) 50, 2695 (1989). 

22. J. A. Holm, Master-Slave Implementation fiir PVM, (Universitat Gottingen re- 
port, 1995, unpublished). |— ■ 

23. Altshuler, Gefen, and Imry claim in Ref.EEl that for the calculation of S^(/i, ip) — 
S^(/i, 0) it is sufficient to consider the function K2{e, e', Lp) only for energies in the 
vicinity of the Fermi energy. Our numerical results show that for non-interacting 
electrons this is certainly not the case. Implicitly these authors seem to assume 
that this property is due to inelastic processes. However, the phenomenological 
inelastic cutoff P that is usually introduced into the Cooperon cannot explain 
this effect. We have not been able to find a microscopic calculation in the litera- 
ture which explicitly shows how dephasing restricts the energies contributing to 
dY?{^, ip)/dp to a small interval close to the Fermi energy. 

24. Ya. M. Blanter, Phys. Rev. B 54, 12807 (1996). 



